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Abstract 



In the stochastic cellular automaton of rule 18 defined by S. Wolfram [Rev. 
Mod. Phys. 55, 601 (1983)] the crossover of the critical behaviour induced by 
nonlocal site exchange was investigated using the coherent anomaly method. 
The continuous variation of the critical (5 exponent was confirmed. 
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The addition of nonlocal exchange to models of short range interaction has recently 
become a major topic of interest. fl| Long-range effective interactions can be defined by 
this kind of dynamics. This is important in simulation because time consuming interaction 
calculations can be substituted by algorithms with simple nonlocal particle exchange. 

We have already shown || that a kind of long-range mixing can be taken into account into 
analytical generalized mean-field (GMF) |||4| calculations as well. The stochastic cellular 
automaton (CA) of rule 18 [BJ with nonlocal site exchange was investigated and the variation 
of the critical point and the order-parameter exponent was determined by extrapolation and 
simple fitting. 

The coherent anomaly method (CAM) || and its enhanced version - introduced very 
recently - have already produced accurate critical exponents for equilibrium statistical 
systems || and for dynamical systems The motivation underlying the present work 

was to show that the GMF + CAM method is able to describe the continuous variation of 
the critical exponent in a dynamical system. The results agree with former simulations and 
GMF extrapolation data. 

Cellular automata models are very useful tools for modelling dynamical systems in many 



branches of science ||14|| . Even one-dimensional models possess very rich features ranging 
from self-organized critical behaviour to chaotic phenomena. Since these models are non- 
equilibrium systems, they can exhibit phase transitions even in one dimension. 



The basic model to be investigated was the one-dimensional, two-state, rule 18 CA [14 
with probabilistic acceptance rate : 

r x if s (t,j - 1) = o, 8 (tj) = o, s(t,j + 1) = i, 

s(t + 1, j) = I X if s(t,j - 1) = 1, s(t,j) = 0, s(t,j + 1) = 0, 
[ otherwise 

where X e {0, 1} is a two valued random variable such that 

Pr(X = l)=p . 

This synchronous update can be followed by a sequential site exchange rule : mx the number 
of living cells ('ones') are selected randomly and swapped with other site values (either zeros 
or ones) selected at random from all over the lattice. This rule preserves the number of 
living cells. 

Taking the steady-state concentration of ones as the order-parameter c, the stochastic 
model without site exchange exhibits a continuous phase transition to an empty state if the 
acceptance probability p is less than p c . This critical probability was found by steady-state 
simulation to be p c = 0.8086(2) 0. Time dependent simulation fTS! - considered to give 



more accurate results owing to the lack of finite size effects - gave p c = 0.8094(2) jTjJ. Both 
simulations predicted that the universality of the critical transition belonged to the directed 



percolation (DP) class flTTI]. The order-parameter exponent of this class is (3 ~ 0.2769(2) 



[0|. Simple mean-field calculation results in p^ IF = 0.5 and (3mf = 1- The generalized 



mean field calculation, which is based on the solution of steady-state equations for n-block 



probabilities, gave a series of approximations This CA in the steady-state can be 



mapped to a simpler rule (rule 6/16) with the new variables 01 — > 1 and 00 — > 0: 
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t-1: 00 01 10 11 

t : 1 10 

and the GMF equations can be set up by means of pair variables. Figure [TJ shows the phase 
diagram of the model calculated by the GMF method and simulation. Pade extrapolation for 
the n = 6 level GMF solution resulted in p c = 0.7986 and (3 = 0.29 |||]. CAM extrapolation 
for the n = 7 pair GMF solution gave an estimate of (3 — 0.2796(2) [p"3| . 

If this stochastic CA model is supplemented with site exchange, the correlations are par- 
tially destroyed and the transition behaviour is shifted towards the simple classical mean-field 
approximation results. We thus see a crossover of the universality. Steady state simulation 
for nearest neighbour site exchange (local) and nonlocal exchange resulted in a similar con- 
tinuous crossover as a function of mixing strength, but on different scales [§]. Though it is 
understandable that the nonlocal exchange is much more effective than the local mixing, it 
is not obvious for a nonlocal site exchange - where a particle can jump to any distance - 
that the crossover is still continuous. To check this, simulations on different sized lattices 
and GMF calculations were carried out. The mapping onto pair variables is not possible 
in this case, since site exchange destroys the ordered structure of the steady-state that - it 
cannot be built up from 00 and 01 pairs. GMF calculations up to the n = 5 level point 
correlation were performed; these were followed by simple extrapolations ||. The variation 
of the p c {iTi) and the (5{m) was found to be continuous in agreement with simulations. In 
what follows I show that an independent extrapolation method (CAM) predicts very similar 
results for the critical exponent, based on finite size scaling . 

Calculating critical exponents of second order phase transitions is a challenging problem. 
Besides simulation there are a number of analytical methods. Theoretical tools for studying 
nonequilibrium statistical physics processes are under development. The coherent anomaly 
method introduced by Suzuki |J is based on a series of generalized mean-field approximations 
and uses finite size scaling theory to estimate the quantities Q n that depend on long-range 
correlations. 

Quantities of the n-th level of approximation in the vicinity of the critical point can be 
described by the classical singular behaviour multiplied by an anomaly factor (a(n)) : 

Q„~a(n)(p/j£-l) w - , (1) 

where p is the control parameter and u c i is the classical critical index. The anomaly factor 
diverges in the n — > oo (and p™ — > p c ) limit, but scales as : 

a(n) ~ (p n c - p c y~^ (2) 

thereby permitting the estimation of the true critical exponent uo. 

Many solved and unsolved equilibrium and nonequlibrium critical systems have been 
studied by this method with success ||. Recently a new parametrization was suggested 
instead of using control parameter p, 

5n = (Pc/P n c ) 1/2 - (P n c /Pc) 1/2 , (3) 

that has an invariance property: p <-» p^ 1 . Taking into account correction terms the anomaly 
scaling now reads as : 
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a(n) = b <^- & + c 5%-^ +1 + 



(4) 



Using this 'enhanced' CAM, very accurate critical indexes were found for the 3 dimen- 
sional Ising model @ . The DP universality of the nonequilibrium stochastic rule 18 CA was 
confirmed |Tj| with accurate f3 exponent estimation. 

Here I used our n = 1...5 point level GMF approximation data with the effect of nonlocal 
site exchange included. As was explained in [[J this mixing with m « 1 can be taken into 
account in the steady-state block probabilities by considering the effects of single particle 
"jump in" or "jump out" of a given n-block. For greater m the mixing can be well described 
by iterating the one-particle jump effect on the equations k = mjnv times such that rrv 
is small. This kind of calculation is in agreement with the sequential simulation process. 
Our experience was that the selection of rrv < 0.001 did not improve the results therefore 
nv = 0.001 was fixed through the calculations. 

The anomaly factor was determined for each m for n = 1...5. 

a(n,m) = c(n,m)/(p/p c (m) — l) /3jv/F , (5) 

where c is the concentration, (3m f = 1- The exponent (3 was fitted out with non-linear 
regression using formula (£|) and the p c (m) values of simulation. Alternatively one could use 
p c {m) estimates determined from extrapolation of GMF data. The stability of the fitting was 
checked by selecting different subsets : n = {1,2,3,4,5}, {1,2,3,4}, {1,2,5} and {3,4,5} 
of the GMF data. The uncertanity was only a few percent for this model. 

In Figure |2| I have plotted the anomaly factors as functions of m. The n — 1 classical 
mean- field approximation factor does not feel the effect of mixing and remains constant. This 
is understandable since in the classical mean-field approximation correlations are neglected 
and the mixing cannot destroy them. Higher levels of approximation for m > 1 give the 
same anomaly factor value, which means according to eq.([|) that the true critical behaviour 
is described by (3 = 0mf = 1- As the mixing decreases from m = 1 the anomaly factors 
of n > 1 change continuously and so does the exponent (3. At m = the anomaly fators 
a(l, 0), a(3, 0), a(5, 0) coincide with the anomaly factors of the pair-correlation GMF results 
a(l), a(2), a(3) |lj. 

Figure |3] shows that the results obtained from this CAM calculation agree very well with 
steady-state simulations and with earlier calculations @. The crossover from [3mf = 1 
begins at m ~ 1. The universality changes continuously toward the DP class in the m — » 
limit (characterised by (3dp — 0.28 ). 

In conclusion, the efficiency of the CAM method for describing universality change in a 
nonequlibrium model is proven by this work. The agreement with other results strengthens 
the belief that nonlocal exchange causes continuous crossover in this model. This method 
shows a possibility for studying models with long-range effective interactions, generated by 
site exchange dynamics. 

This research was partially supported by the Hungarian National Research Fund (OTKA) 
under grant numbers T-4012 and F-7240. 
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FIGURES 



FIG. 1. Phase diagram of the stochastic rule 18 CA determined from simulation (diamonds) 
and GMF pair- approximations at levels n = 1, 8 

FIG. 2. Anomaly coefficients as functions of mixing strength m obtained from 
n= 1(A), 2(x), 3(a), 4(+), 5(o) levels of GMF approximations. 

FIG. 3. Crossover induced by nonlocal site exchange m from DP to mean-field universality. 
The phase diagram was calculated by CAM (+) and steady-state simulation (x). The DP value 
of p c ~ 0.276 is also shown. 
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